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We present a numerical study of conductance oscillations of transition metal multilayers as a func- 
tion of layer thickness. Using a material-specific tight-binding model, we show that for disorder-free 
layers with random thicknesses but clean interfaces, long-period oscillations in the conductance can 
occur, which are reminiscent of those found in structures exhibiting GMR. Using a heuristic effective 
mass model, we argue that these oscillations arise from beating between the Fermi wavevector and 
a class of wavevectors characteristic of the superlattice structure. 



PACS: 73.23.Ad, 73.40.Jn, 73.61.At 

Oscillations in transport properties of metallic su- 
perlattices have been largely studied in magnetic/non- 
magnetic multilayers, which exhibit giant magnetoresis- 
tance (GMR) ||-||]. In such systems the magnetic config- 
uration of the superlattices is an oscillating function of 
the layer thickness, switching periodically between fer- 
romagnetic and antiferromagnetic alignment of adjacent 
magnetic layers, with typical periods extending over sev- 
eral atomic planes. This arises from an oscillating inter- 
layer exchange coupling and results in oscillations of the 
overall conductance, with the antiferromagnetic configu- 
ration being the more resistive. 

Recently a new set of measurements on Ni/Co 
Ag/Pd 0, Ag/Au and Ag/Cu Hi multilayers revealed 
the possibility of long-period oscillations of a different 
origin. On the one hand, the Ag based multilayers are 
entirely non-magnetic. On the other, the Ni/Co multi- 
layers were measured in high magnetic field, far above the 
coercive field of the structure, which rules out magnetic 
misalignment between magnetic layers as source of the 
oscillations. In these experiments, all the measurements 
were conducted with the current in plane (CIP) config- 
uration and the authors associated the oscillations with 
the formation of a d superlattice bound state, giving rise 
to a strong s-d scattering. 

In this letter we predict that such oscillating behaviour 
can also occur with the current perpendicular to the 
planes (CPP), in clean superlattices with very good inter- 
faces, but where the layer thickness fluctuates randomly. 
To the best of our knowledge this is the first time that 
conductance oscillations in the CPP configuration have 
been identified. 

To address this problem, we have developed a very 



efficient technique to calculate transport properties of a 
finite multilayer attached to semi-infinite pure crystalline 
leads, as sketched in figure |l|. Our calculations are based 
on the Landauer-Biittiker formalism using ab initio 
s-p-d tight-binding Hamiltonians with nearest-neighbor 
hopping. The conductance F"' of a given spin species is 
obtained from the Landauer formula M 



h 



(1) 



where T*^ is the total transmission coefficient for the spin 
a ((T =ti i) calculated at the Fermi energy. The latter is 
obtained by extracting the S matrix from the total Green 
function G of the superlattice in contact with external 
leads. The total Green function is calculated via Dyson's 
equation, starting from the surface Green function g of 
the leads and an effective Hamiltonian Hcs describing 
the finite multilayer. A detailed description of the tech- 
nique will be given elsewhere joj. In all the following 
calculations we consider a perfect lattice match between 
clean fee layers and hence fc|| is a good quantum number 
(the symbol || represents the in-plane coordinates and 
the symbol _L the direction of the current perpendicular 
to the planes). The Hamiltonian is diagonalised in the 
Bloch basis to yield 



(2) 



where the sum over fc|| extends over the two dimensional 
Brillouin zone. In what follows, we employ of order 10^ 
/c|| -points, which is sufficient to render effects due to the 
finite number of fcy-points negligible compared with the 
oscillations of interest. In what follows, for Ni/Co and 
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Ag/Pd multilayers, we calculate the total conductance 
of the two independent spin channels as a function of 
layer thickness, in the limit that the spin-flip and phase- 
breaking lengths are infinite. It should be noted that the 
majority bands of Ni and Co arc s-p-likc and are closely 
aligned. On the other hand the minority bands are d- 
like and possess a relative shift in energy of about 0.7 
eV. Hence we expect a large contribution to the conduc- 
tance from the majority channel and a small contribution 
from the minority channel. For Ag/Pd the situation is 
qualitatively different, because at the Fermi energy the 
DOS of Ag is dominated by s-p electrons, while in Pd it is 
dominated by d electrons. As a consequence one expects 
strong interband scattering at the interfaces between the 
different metals. 

Following reference |]l0|| , we consider a pseudorandom 
layer arrangement, in which a finite A/B multilayer, at- 
tached to semi-infinite leads of material A, possesses B- 
layers of fixed thickness Ib and A-layers of random thick- 
nesses I A which are allowed to fluctuate by ±1 atomic 
planes (AP) around a mean value Ij^. In all the follow- 
ing simulations, we consider multilayers consisting of 10 
A/B bilayers and for each Ib, show results for the average 
conductance of 10 random configurations of the A-layers. 

For Ni/Co and Ag/Pd respectively, figures || and ^ 
show the mean conductance as a function of Ib, along 
with the corresponding root-mean-square deviation from 
the mean. These figures demonstrate the presence of 
long-period oscillations, with amplitudes not exceeding 
25% of the mean conductance. Moreover the Ni/Co sys- 
tem shows smaller oscillations than the Ag/Pd system, 
and despite the fact that the conductance of the majority 
spin channel is almost double that of the minority, the 
oscillations arise predominantly from the minority spins, 
where the scattering is strongest. Bearing in mind that 
these calculations involve the CPP configuration, the size 
of the oscillations compared with the mean conductance 
for the Ni/Co system, is consistent with some of the ex- 
perimental data ^. 

To understand how quantum interference of the con- 
duction electron wave-functions can give rise to such long 
period oscillations, we now develop a heuristic continuum 
model, within the effective mass approximation, describ- 
ing an infinite 3D superlattice with a Kronig-Penney po- 
tential and a parabolic band. The spin-dependent Hamil- 
tonian for such a system is 
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m*{z) dzm*(z)dz 



+ V''{z), (3) 



where Vj„ is the 2D Laplacian. Since the structure of 
possesses translational invariance in the x-y direc- 
tions, the spin-dependent Kronig-Penney potential V^lz) 
and the effective mass m*{z) are functions of z only. 
Consequently the problem can be mapped onto a 
dependent ID problem, whose Hamiltonian is 



H'^{z-k 



d 



2 1,2 



2 diZm*{z)Az 2m* {z) 



+ V''{z). (4) 



For each fc|| and spin tr, an eigenstate at the Fermi energy 
contributes jh to the conductance of this infinite peri- 
odic structure. In the general case, the eigenstates can 
be obtained numerically using standard transfer matrix 
techniques. First consider the case of constant m*(z), 
where the problem can be solved analytically. Since the 
Hamiltonian depends on fc|| only through an energy 
shift, one finds that the conductance per unit area has 
the simple form 



(5) 



where A„ is the bandwidth of the n-th energy band of 
the Hamiltonian 



H%z) = - 
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2m* dz2 



+ V-{z), 



(6) 



and the sum is over all occupied minibands. Consider an 
infinite superlattice composed of materials A and B, with 
layer-thicknesses I a and Ib {I a + Ib = L), and Kronig- 
Penney potential V = Vo {Ep > Vo) in the metal A and 
V — in the metal B. If k± is the Bloch vector in the 
direction of the current, the secular equation is 



cos{k±L) = cos(A;a^a + ^b^b) — 

(/ca + ^b)^ . , X . /, , ^ 
; — ; sm(fcAiA) sm(A;BiB), 

fcAKB 



(7) 



with kA{E) = ^2m*{E -Vo)/h and kB{E) = 
V2m* E/h. Based on this expression, we now argue that 
the bandwidths exhibit several periods of oscillation as 
the layer thicknesses are varied. 

To describe Ni/Co multilayers, we vary the thickness 
of metal B keeping fixed the thickness of metal A. To un- 
derstand the oscillatory behaviour of the band-widths, 
we note that equation (|^) cannot be satisfied at energies 
for which 



kA{E)lA + kB{E)lB — mn, 



(8) 



where m is an integer. Hence at E = Ep and fixed Ia, 
successive bandgaps appear at the Fermi energy Ep when 
changes by 
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Equation introduces the first period of oscillation l^. 
The second period corresponds to the presence of narrow 
gaps below the Fermi energy. From equations (Q) and (||) 
narrow bandgaps appear at the energies 



2m*l\ 



(10) 
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whenever the lengths Ib equal 



/(") _ 



2m*E 



A 
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The total bandwidth A and hence the conductance per 
unit of area (H) are oscillating functions with periods Zg 

and the Zb"'''s. All these periods are of order \p (ie few 
A) , but beating between them can give rise to long-period 
oscillations. It is important to note that the Fermi pe- 
riod is defined only through the Fermi energy, while the 
periods Zg""* depend critically on the superlattice geom- 
etry. In particular, because the energies corresponding 
to periods (^0|) depend on and must not exceed the 

(n) 

Fermi energy, the number of Zg s depends on the thick- 
ness of the metal A. If Za is large, a large number of Zg"^ 
periods will be present and the beating pattern will be 

in) 

complex. On the other hand, if I a is small, few Zg 's will 
be present, giving rise to a simple beating pattern. A nu- 
merical evaluation of Eq. ^ is shown in figure ^. For the 
chosen parameter in this plot, we expect only one Zg"^ and 
clear beats are observed, with period 2Zg^Zg/(/g^ — Zg). 

Since the Zg"'' periods are characteristic of the superlattice 
structure we predict that the period of the long oscilla- 
tions can be set by choosing the appropriate superlattice 
geometry. 

The above dependence of oscillations on the multi- 
layer structure is missed by a trilayer quantum well ap- 
proach to conductance oscillations and GMR , where 
only two periods have been identified. The first of these 
p^^ depends on the extremal Fermi surface radius of the 
spacer forming the well, and in the parabolic band ap- 
proximation corresponds exactly to the period Zg. The 
second period p'^P depends on the cut-off of the sum over 
the fc|| 's, and in the parabolic approximation, on the en- 
ergy difference between the Fermi energy and the step 
potential Vo- In our superlattice description, this period 
is replaced by the class of periods Zg"^ , which are a func- 
tion of the superlattice structure itself. This structural 
dependence of the oscillation periods is the key to under- 
standing the apparent non-reproducibility of the long pe- 
riod oscillations from sample to sample, observed in some 
of the experiments ||^ . It may be shown that these beat- 
ing features are preserved when a more realistic material- 
dependent effective mass is used. Futhermore, for those 
cases where we expect the mass to be significantly differ- 
ent in the two materials (e.g. mainly s-p like in one ma- 
terial and mainly d-like in the other) we have shown that 
the Kronig-Penney model reproduces the main features 
of the more accurate tight-binding model, with physi- 
cally reasonable choices of the band offsets and effective 
masses. 

Bearing in mind that our analysis describes the CPP 
configuration, we can also speculate on the absence of the 
oscillations in other recent experiments Ag/Cu ||] 



exhibits very good phase separation between the differ- 
ent metals and hence it should be a good candidate for 
observing conductance oscillations. However the band 
match between Ag and Cu is very good, resulting in a 
very small scattering potential at the interface. In the 
effective mass approach this means a very small step po- 
tential Vo with respect to the Fermi energy. A large num- 
ber of periods Zg will be present and the beats will be 
difficult to detect. The same argument is valid for Ag/Au 
. In addition the high miscibility of Ag and Au results 
in dirty interfaces. Ag/Pd is in theory a good candi- 
date to show conductance oscillations because of the large 
mismatch between the Ag and Pd bands. Unfortunately 
interdiffusion at the interface is difficult to avoid and the 
elastic mean free path will be quite short. Finally, we 
observe that for Ni/Co the majority band repro- 

duces roughly the situation of Ag/Cu, while the scat- 
tering in the minority band is quite large. According 
to the effective mass model the minority band will pos- 
sess a low conductance with large oscillations, while the 
conductance of the majority band will be large and the 
oscillations small. This is precisely what we obtain from 
in the material specific tight-binding calculations. 

In summary, we have investigated the possibility of 
large long-period oscillations in metallic superlattices in 
the ballistic regime. With the current perpendicular 
to the plane and superlattices with pseudorandom layer 
thicknesses, the oscillations are predicted by accurate 
tight-binding calculations. An effective mass analysis 
provides a qualitative understanding of the nature of the 
oscillations and highlights their dependence on the su- 
perlattice geometry. 
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FIG. 1. Finite multilayer connected to pure crystalline 
seminfinite leads, g are the surface Green function describing 
the leads and Hes is the effective hamiltonian describing the 
multilayer. 
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FIG. 3. Conductance of Ag/Pd multilayers as a function 
of the Pd thickness with an average Ag thickness of 5 atomic 
planes. 
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FIG. 2. Conductance of Ni/Co multilayers as a function of 
the Ni thickness. The Co thickness is 10 atomic planes. The 
inset shows the two spin conductances on the same scale with 
the upper plot for majority spin and the lower for minority 
spin. 
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FIG. 4. Conductivity per channel in the effective mass 
approximation. The parameters are Ef = lOeV, Vo = 6eV, 
m* = 0.5MeV, Ia = SA. 
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